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Abstract 

In this paper, the physical realizability property is 
investigated for a class of nonlinear quantum systems. 
This property determines whether a given set of non- 
linear quantum stochastic differential equations corre- 
sponds to a physical nonlinear quantum system satisfy- 
ing the laws of quantum mechanics. 



1. Introduction 

In the physics Hterature, methods have been devel- 
oped to model a wide range of open quantum systems, 
such as those encountered in quantum optics, within 
the framework of quantum stochastic differential equa- 
tions O, (gl, lEl. In this context, quantum noise is used 
to represent the influence of large heat baths and bo- 
son fields, including optical and phonon fields, from 
which completely positive maps, Lindblad generators 
and master equations are obtained by taking expecta- 
tions ifTOl . We may then distinguish two types of quan- 
tum stochastic differential equations, namely, linear and 
nonlinear. Linear quantum systems can be described 
by linear quantum stochastic differential equations and 
arise mostly in the area of quantum optics; e.g., see 
Q, ||9]| and Q. An important class of linear quantum 
stochastic models describe the Heisenberg evolution of 
the annihilation and creation operators of several in- 
dependent open quantum harmonic oscillators that are 
coupled to external coherent bosonic fields, such as co- 
herent laser beams; e.g., see O, Q and ||9l. A spe- 
cial class of these linear quantum systems is driven by 
quantum Wiener processes as in O where physical re- 
alizability conditions are developed to determine when 
the linear quantum system under consideration can be 
regarded as a representation of a linear quantum har- 
monic oscillator. This imposes some restrictions on the 
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matrices describing the linear quantum model. This no- 
tion of physical realizability has been further investi- 
gated in 1 11 1 where the authors prove the equivalence 
between the algebraic conditions for physical realiz- 
ability obtained in | 2| and a frequency domain condi- 
tion that an associated linear system is (/,/) -unitary. In 
fact, fTT] extends the frequency domain physical realiz- 
ability results of lfT2l . llT3l and 1 1 1 which apply to linear 
quantum systems described purely in terms of annihi- 
lation operators only to a more general class of linear 
quantum systems, which are described in terms of both 
the annihilation and creation operators. More explic- 
itly, in lfT2l . lfT3l and O, the physical realizability con- 
ditions developed for annihilation-operator linear quan- 
tum systems have been related to the lossless property 
of linear systems and |11| generalizes this result to re- 
late the corresponding physical realizability conditions 
of annihilation and creation operators linear quantum 
systems to the property of (/,/) -unitary. However, the 
question of addressing the physical realizability condi- 
tions for nonlinear quantum systems described by non- 
linear quantum stochastic differential equations is still 
open. And in this paper, we restrict attention to a class 
of nonlinear quantum systems which is a generalization 
of the annihilation operator only linear quantum sys- 
tems considered in 11121 . 11131 and Q. 

The paper is organized as follows: Section II de- 
fines the class of nonlinear quantum systems under con- 
sideration. Section III develops conditions under which 
the commutation relations for the nonlinear quantum 
systems are preserved. Section IV defines physical real- 
izability in terms of nonlinear open quantum harmonic 
oscillators. Section V provides an example to illustrate 
the theory developed. Section VI concludes the paper. 

2. The Class of Nonlinear Quantum Sys- 
tems 

We consider an open quantum system G with phys- 
ical variable space sz/q consisting of operators defined 



on an underlying Hilbert space M'g- The self-energy 
of this system is described by a Hamiltonian G s^q- 
The system is driven by m field channels given by the 
quantum stochastic process W 



These describe the annihilation of photons in the field 
channels and are operators on a Hilbert space F , with 
associated variable space ^ . In that case, F is the 
Hilbert space defining an indefinite number of quanta 
(called a Fock space | 5 1), and ^ is the space of opera- 
tors over this space. 

We assume the process W is canonical, meaning 
that we have the following second order Ito products: 



dw^(t)dmi{tY 
dw^(tYdWi(tr 



Skidt; 
0; 
0; 




where Wk{ty is the operator adjoint of Wk{t) which is 
defined on the same Fock space. 

The system is coupled to the field through a scat- 
tering matrix S = 1 and a coupling vector of operators L 
given by 



V L„ 



J 



where Lj G s^g- 

The notation G = {S^L^H) is used to indicate an 
open quantum system specified by the parameters S^L 
and H where H represents the Hamiltonian of the sys- 
tem. For S = I, the Schrodinger equation is given by 

dU{t) = ^^dW'^L-LUw - ]^L^Ldt - iHdt^U(t) 

(1) 

with the initial condition U (0) = /. Equation ([T]) deter- 
mines the unitary motion of the system in accordance of 
the fundamental laws of quantum mechanics. Note that 
the notation refers to the Hilbert space adjoint. 

Given a system annihilation operator ai G s^g^ its 
Heisenberg evolution is defined by ai{t) = jt{ai) = 
U{tyaiU{t) and satisfies 

dai{t) = {^L{t){ai{t))-i[ai{t),H(t)]) (2) 
^[L(t)\ai(t)\dW(t) 

where the notation [A^B] = AB — BA denotes the com- 
mutator of two operators A and B. 



In equation dJ]), all operators evolve unitarily and 
the notation refers to 

^^(a/) = V[^/,L] + i[L^a/]L. (3) 

In this paper, we restrict attention to annihilation only 
coupling operators; i.e, the coupling operator L(t) de- 
pends only on the annihilation operator ai{t) and not 
the creation operator a/ (^)*. Therefore, L{t) satisfies the 
following commutation property: [a/,L] =0 and then 

^L{ai)=^-[L\ai\L. (4) 

The generator of the system G is given by 

^G{ai) = -i[ai,H]^^L{ai), (5) 

For the case of having n annihilation operators 
. . . we define 

/ ai 



\ an 



^LM) = \[L\ai]L 



(6) 



and 



^G,(^) = -/[^,-,if]+^L,(a). (7) 

Therefore, we can write 

dai{t) = {^L,(t){ai{t))-i[ai(t),H{t)]) 
^[L{t)\ai(t)\dW(t) 

and 

da{t) = {^L{t){a{t))-i[a{t),H{t)]) 
-^[L{t)\a{t)]dW{t) 



where 



and 



[a,H] 




(8) 



(9) 



Let 



A{a{t),a{ty) 



I Ai{a{t),a{t)^) 
\ An{a{t),a{ty) 
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/ [auH\ 
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\ [an,H] 


J 



^L{a)-i[a,H] 



and 



B{a{t),a{t)^) 



I Bi{a{t),a{ty) 
\ Bn{a{t),a{ty) 



(10) 



(11) 



= [L\a]. (12) 

Hence, 

da{t) = A{a{t),a{t)'^)dt^B{a{t),a{t)'^)dW{t) 
da(ty = A{a{t),a{tYydt^B{a{t),a{tyydW{ty. 

Note that for the case of matrices, the notations * and 
refer respectively to the complex conjugate and the 
complex conjugate transpose of the matrix in question. 

(13) 

On the other hand, the components of the output fields 
are defined by y{t) = jt{W{t)) = U(tyw(t)U(t) and 
satisfy the nonlinear quantum stochastic differential 
equations 

dy{t) = C{a{t))dt^D{t)dW{t) 

dy{ty = c{a{t)ydt^D{tydw{ty 

(14) 

where C{a{t)) = L{t), C{a{t)y = L{ty, D{t) = I and 
D{ty = I. Hence, the system G can be described by 
the following nonlinear quantum stochastic differential 
equations 

da{t) = A{a{t),a{ty)dt^B{a{t),a{ty)dW{t)\ 
dy{t) = C{a{t),a{ty)dt^D{a{t),a{ty)dW{ij5) 



where a{t) = 
C{a^a^) = 



a 

a* 



A{a^a^) 



C{a) 
C{ay 



A{a^a^) 
A(a^a^) 

B{a,a^) = 



dW{t 



B{a,a^) 
B{a,a^) 
dW{t) 

dw{ty 



,D{t): 



D{t) 
D{ty 



and 



Definition 2.1 The class of nonlinear quantum system 
we consider in this paper are nonlinear quantum sys- 
tems that can be represented by the QSDES f [75l ) such 
that the matrices A{a,a'^) and C {a) satisfy [C{a)^a^] = 
0, [A(a,a'^),a^] = — [a,A(a,a'^)^] and 



Ai{a,a'^) 
Cy{a) 



E E £ £ apik,h/-{a*i f' 

ki=Ohi=Ol=lp=l 
/t„=Op=l 



(16) 



for i = 1 , . . . , /t and v = 1 , . . . , m with m^., m^., m^^, fJ^^p 
and CCpikihi integers. 

Moreover, the matrices A{a^a^), B{a^a^) andC{a) 
satisfy the following property: 



1 

2n 



A{a,a^)^ 6 ^a^a 



In 



- [a^Q ^A{a,a^),a\ 



A(a^a^) — -B{a,a^)C{a,a^) 



(17) 



where n = sup {ki + hi) 

\<i<n 





0* 



and = For more information, on 0, refer to 

Section III below. 

In the sequel, we assume that dW{t) admits the fol- 
lowing decomposition: 

dW{t) = P,y{t)dt ^ dw{t) 

where w{t) is the noise part of W{t) and I5w{t) is an 
adapted process (see f5\ and f6l). The noise w{t) is a 
vector of quantum Weiner processes with Ito table 



dw{t)dw{t) 



t . 



■ F^dt 



(18) 



(see [5]) where Fy^ is a non-negative definite Hermitian 
matrix. Here, the notation represents the adjoint trans- 
pose of a vector of operators. Also, we assume the fol- 
lowing commutation relations hold for the noise com- 
ponents: 



[dw{t),dw{ty] =dw{t)dw{ty -{dw{tydw{tff ■■ 

(19) 

Here 7^ is a Hermitian matrix and the notation ^ de- 
notes the transpose of a vector or matrix of operators. 
The noise processes can be represented as operators on 
an appropriate Fock space; for more details, see Q. 



■- T^dt. 



The process represents variables of other sys- 
tems which may be passed to the system ([TSl) via an 
interaction. Therefore, it is required that Pw{0) be an 
operator on a Hilbert space distinct from that of gq and 
the noise processes. We also assume that Pw{t) com- 
mutes with a{t) for all ^ > 0. Moreover, since liw{t) is 
an adapted process, we note that Pw{t) also commutes 
with dw{t) for all ^ > 0. We have also that [a, dW^] = 0. 

3. Commutation Relations 

For the nonlinear quantum system ([15]), the initial 
system variables a{0) = ao consist of operators satisfy- 
ing the commutation relations 

[aj{0),al{0)] = Qjk, j,k=l,...,n. 
[aj{0),al{0)] = ejk, j,/:=l,...,2^.(20) 

Here, the commutator is defined by — cijal — 

d^aj = (djk where is a complex matrix with elements 
(djk. With — {ai^...,an), the relations (l2Qb can be 
written as 



(2, a' 



aa 



\a,a 



t1 A 





0* 







= 0. 



(21) 



3.1. Preservation of the Commutation Rela- 
tions 

The following theorem provides an algebraic char- 
acterization of when the nonlinear quantum system 
([T5l) preserves the commutation relations (l2Qb as time 
evolves. 

Theorem 3.1 For the nonlinear quantum system 
( [73T) . we have that [ap(0),^*(0)] = 0^^ implies 
[dp{t)^dq{t)] = for all t > with p^q = 1 . . . ,2^ 
if and only if 

\_A{d,d^),d^\ + \d,A{d,d^)'^\ 
-^B{d,d^)T^B{d,d^)^ =0; 
[B{d,d^),d^] = and 
[d,B{d,d'^)'^] =0 



(22) 



with Tw 











4. Physical Realizability And The Nonlin- 
ear Quantum Harmonic Oscillator 



sponds to a nonlinear open quantum harmonic oscilla- 
tor. Such a system is said to be physically realizable. 
The class of nonlinear open quantum harmonic oscil- 
lators under consideration are defined by the nonlinear 
coupling operator L and a nonlinear Hamiltonian H. To 
derive a nonlinear system of the form ilSh from a non- 
linear open quantum harmonic oscillator defined by L 
and H, we proceed as in Section II to get the following 
definition. 

Definition 4.1 A nonlinear quantum system of the form 
f/TI) is said to be physically realizable if it is a repre- 
sentation of a nonlinear open quantum harmonic oscil- 
lator defined by a nonlinear coupling operator L and 
a nonlinear Hamiltonian H, i.e., if there exist H and 
L such that the matrices A{a^a^)^B(a^a^)^C{a^a^) and 
Dia^a^) satisfy the following equations 



A{a^a^) 

B{a,a^) 
C{a,a^) 
D{a^a'^) 



i [L^ ,d]JL^i[H,d]; 



[L\d\ /; 

u 



(23) 



where J - 



I 

-/ 



The following theorem provides necessary and suffi- 
cient conditions for a nonlinear quantum system of the 
form ([T5l) to be physically realizable. 

Theorem 4.1 A nonlinear quantum system of the form 
f/5l) is physically realizable with T^^ = J - 
if and only if 







[A(a,a'^),a 



^ [d.Ai^a.a'^)^] 
^B{a,a^)jB{a,a^)^ 



[B{a,a^),d^] =0; 
[d,B{a,a^)^] =0; 
B{a^a^) — \C{a.,a^)\d\j 
D{a,a^)=I. 



and 



(24) 



In this case, the corresponding nonlinear self-adjoint 
Hamiltonian H is given by 



H= — {d'^0-^A{a,a'^) -A(a,a'^)'^e"^a) (25) 



In this section, we consider conditions under which 
a nonlinear quantum system of the form STSh corre- 



where n is as defined in Definition \2.1\ and the corre- 
sponding coupling operator is L = C{a^a^). 




Figure 1. Schematic for tlie control of an opti- 
cal parametric oscillator. 

5. Example 

We consider the application of the physical real- 
izability conditions developed in this paper to an opti- 
cal squeezer. A pulsed laser is used to drive the optical 
squeezer. The use of such a pulsed laser provides many 
advantages in an experimental situation; e.g, see |4|. 

The squeezer under consideration is implemented 
as an OPO driven by two time- varying optical fields; see 
Figure 1 . A full description of the apparatus of Figure 1 
is described in lfT4l . 

5.1. Modeling 

The dynamics of the OPO system can be described 
as in 041 but with time-varying matrices and fields 
reflecting the fact that a pulsed laser is used to drive 
the experiment. The dynamics of the optical subsys- 
tem contain nonlinear terms. The quantities Ain and Bin 
are the input fields which set-up the fundamental and 
second-harmonic intracavity fields ai and a2. Ka-^ and 
Kij.^ are the loss rates of the input/output mirrors for the 
ai and a2 fields respectively. The parameters Ka^ and 
Kij^ are the internal loss rates for the corresponding two 
fields. Also, Ka^ = Ka-^ + and Ka^ = Ki,.^ + K*^^ are 
the associated total resonator decay rates. Furthermore, 
5Ai and 5^/ represent the vacuum fields due to the in- 
ternal losses. The output fields are given by Aout and 

^out- 

The Heisenberg equations of motion are given by 
(see (H): 



da\{t) 
da2(t) 



Ai{a,a^){t)dt ^ Bi(t)dwi(t) 
A2 {a,a^)(t)dt + B2 (t)dw2 {t). 



In that case, A\{a^a^) = —Ka^a\ — 2x^^^a\a2, B\ = 



dwi{t) = —dAin and Jw2(0 = —dBin. a\ anda2 satisfy 
the following commutation relations: 



[ai,a\] 
[a2,al] 

[a2,a\] 
[a\,a2] 



Let a = 





a\ 






CI2 






B2 


, dw 




,«t) 1 

at)* _ 


,B 


' B 




dwi 

dW2 



, a 



Ai{a^a^) 
A2{a^a^) 
a 
a' 



B = 



, A{a^a'^) 



and dw -- 



dw 

dw"" 



On the other hand, we define 

dy{t) = C{a){t)dt^D{t)dw{t) 



where C{a 

Also, we define y 
and 5(0 : 



y 
y 



Sind D{t) 



I 
/ 

C{a) 
C(a)* 



D{t) 
D{ty 

Then , the OPO system can be described by 



da(t) = A{a,a^)(t)dt^B(t)dw(t)\ 
dy(t) = C{a,a^)(t)dt^D(t)dw(t). 

Note that this OPO system fits within the class of non- 
linear quantum systems considered in this paper. In fact, 
the matrices A{a^a^) and C{a) satisfy [C(a),(2^] = 0, 
\A{a^a^)^a^\ = — [a,A(a,a'^)^] and 



1 

2h 



1 



In 



A{a^a^) — -B{a^a^)C{a^a^) 



where ^ = 3,0 = 



In fact, ^ 



/ 

-/ 



and = = /. 



A((2,(2') 6 a, a 



and ^ ^A(a,a'^),a] 



X^^^ala2 



-X^^^cila2 
-X^^^alai 



. On the 



other hand, A (a, a^) — ^B{a^a^)C{a^a^) 



. Hence, 



1 

2n 



In 



In addition, the system under consideration is physi- 
cally realizable. In fact, B{a,a^) = [C{a,a'^y ,a]j = 













-^y: 



^ -^/2k^2 

On the other hand, 

[d,B{a,a'^)'^] = and = /. 

Also, [A{a,a'^),d'^] 

2x^^^ai -ka, 

-Ix'^^^al 

d^Aia^a^) 



0; 



2;t;(2)a* -Ix^^W 


2x^M 






-2x^^^ai 







-2k 



'CL2 



2x^^^a\ 


and 

B{a,a^)JB[a,a^)^ 
2ka, 

2ka^ 
-2ka, 


Hence, [A(a,a'^),a'^] + [(2,A(a,a'^)'^] + 

B{a,a^)JB{a,a^y 

In this case, the corresponding nonlinear self- 
adjoint Hamiltonian H is given by 

H = ^{d'^e-^A{a,a'^)-A{a,a'^Ye-^d) 
= ix^'^afa2-iX^Mal 
and the corresponding coupling operator is L = 

6. Conclusion 

In this paper, the issue of physical realizability for 
a class of nonlinear quantum systems is developed and 
verified by applying the theory developed to an optical 
squeezer. A future extension of this paper would be to 



relate the physical realizability properties to the loss- 
less properties of the class of nonlinear quantum sys- 
tems considered in this paper. 
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